#!/usr/bin/env python

import matplotlib.pyplot as pl
import numpy as np
#data = np.loadtxt('d-distribution-data/d-q6.txt')
data,xxx = np.loadtxt('d.txt', usecols=(1,0), unpack=True)
(hist,bins) = np.histogram(data, bins=100)
#own norming
mysum=float(sum(hist))
hist = hist/mysum
pl.plot(bins[1:],hist)
pl.legend(loc='upper left')
pl.xticks(np.arange(0,1.1,0.2))
pl.yticks(np.arange(0.02,0.1,0.02))
#pl.xlim(0,1.1)
pl.axvline(0.7, lw=1)
#overrideTitle=r'Q6, ts=0.7'
#pl.ylabel(r'P($q_6$)')
#pl.xlabel(r'$q_6$')
pl.show()

